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Supermassive Black Holes in Elliptical Galaxies: Switching 
from Very Bright to Very Dim 
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ABSTRACT 

Relativistic outflows (mainly observed in the radio) are a characteristic feature of both 
Galactic stellar mass black holes and supermassive black holes (SMBHs) . Simultaneous 
radio and X-ray observations of Galactic sources have shown that the outflow is strong 
at low accretion rates, but it weakens dramatically or disappears completely at high 
accretion rates, manifesting structural changes in the accretion flow. It is reasonable 
to assume that SMBHs follow the same trend. For low luminosity SMBHs in nearby 
elliptical galaxies and clusters, recent observations strongly suggest that the outflows 
play the central role in keeping the gas hot (mechanical feedback). If the outflow 
is quenched in SMBHs at high accretion rates similarly to the behavior of galactic 
sources, then the straightforward consequence is a relatively weak feedback of rapidly 
accreting SMBHs. We argue that elliptical galaxies and their central engines should 
then evolve through two stages. Early on, the central SMBH rapidly grows by accreting 
cooling gas at a near-Eddington rate with high radiative efficiency but with weak 
feedback on the infalling gas. This stage terminates when the black hole has grown 
to a sufficiently large mass that its feedback (radiative and/or mechanical), despite 
the low gas heating efficiency, is able to suppress gas cooling. After that the system 
switches to a stable state corresponding to passively evolving ellipticals, when the 
accretion rate and radiative efficiency are very low, but the gas heating efficiency is 
high and energy input from the relativistic outflow keeps the gas hot. 

Key words: galaxies: active - galaxies: jets - galaxies: nuclei 



1 INTRODUCTION 

It is now widely accepted that virtually every elliptical 
galaxy hosts a supermassive black hole (SMBH) at its cen- 
ter with its mass tightly correlated with the mass and stel- 
lar velocity dispersion of the galaxy itself (e.g. Ferrarese 
& Merritt 2000; Tremaine et al. 2002). During the early 
stages of galaxy evolution, when large quantities of cold gas 
are presumably present, these black holes accrete matter 
at high rates and are observed as bright QSOs once they 
have grown to a considerable mass. In contrast the central 
SMBHs of local elliptical galaxies usually have very small 
bolometric luminosities. Recent X-ray and radio observa- 
tions suggest that, despite low observed luminosities, these 
SMBHs in local ellipticals eject material outflows that have 
a strong impact on the thermal state of the ambient hot 
gas, which is characterized by a relatively short cooling time 
and should cool and condense in the absence of an external 
source of energy. These outflows usually manifest themselves 



as bubbles of radio emitting plasma coincident with depres- 
sions in the X-ray images (e.g. Bohringer et al. 1993; Fabian 
et al. 2000; McNamara et al. 2000). The energetics of the 
outflows/bubbles and their evolution have been extensively 
discussed (e.g. Owen, Eilek & Kassim 2000; Churazov et 
al. 2001, 2002; Briiggen & Kaiser 2001; David et al. 2001; 
Nulsen et al. 2002; Fabian et al. 2003; Ruszkowski, Briiggen 
& Begelman 2004; Birzan et al. 2004). At least in several 
well studied objects the power of the outflow exceeds the 
observed bolometric luminosity and is sufficient to offset the 
gas cooling losses. 

Radio emitting flows are also observed in X-ray bina- 
ries containing stellar mass black holes and the physics there 
is believed to be similar to that in SMBHs. For individual 
binaries the radio emission has been observed to vary de- 
pending on the X-ray flux and on the spectral state of the 
source (e.g. Fender & Belloni 2004). In particular it was 
found that an outflow ceases above a certain accretion rate 
(e.g. Fender et al. 1999). There are indications that active 



© 2001 RAS 



2 Churazov, Sazonov, Sunyaev, Forman, Jones and Bohringer 



galactic nuclei (AGNs) may also exhibit a similar kind of 
radio power behavior (Maccarone, Gallo & Fender 2003). It 
thus seems reasonable to assume that the outflow behavior 
is similar for galactic binaries and SMBHs and consider the 
most straightforward implications for elliptical galaxies and 
their central SMBHs. This is the purpose of the present pa- 
per. Previously the role of feedback from growing SMBHs on 
their host galaxies has been discussed by many authors us- 
ing different assumptions (e.g. Binney & Tabor 1995; Silk & 
Rees 1998; Ciotti & Ostriker 2001; Cavaliere, Lapi & Menci 
2002; Wyithe & Loeb 2003; Granato et al. 2004; Murray, 
Quataert & Thompson 2005; Sazonov et al. 2005; Springel, 
Di Matteo & Hernquist 2005, Begelman & Nath, 2005). 



2 BLACK HOLE ENERGY RELEASE IN 
X-RAY BINARIES AND AGNS 

Simultaneous observations of several X-ray binaries, in par- 
ticular black hole candidates, in radio and X-ray bands have 
recently led to the following broad picture (e.g. Gallo, Fender 
and Pooley 2003; Fender & Belloni 2004). In the so-called 
"low/hard" source state X-ray and radio fluxes are well cor- 
related with a functional form F r oc F x ' 7 . In this state the 
accretion rate is low - M <C M E dd, the accretion flow is 
optically thin (geometrically thick) and radiative efficiency 
(luminosity/Mc 2 ) is much lower than canonical 10%. Here 
M E dd = 47rG ^T <fgH ' ^ BH i s tne black hole mass, G gravita- 
tional constant, m p proton mass, c speed of light, or Thom- 
son cross section, and r\ ^10%. For individual sources (a no- 
table example is GX339-4) the F r /Fx correlation extends 
over at least three orders of magnitude in X-ray luminosity. 
It can be extended by another 7-9 orders of magnitude by 
adding a large sample of individual AGNs (Maccarone et al 
2003; Merloni, Heinz & Di Matteo 2003). In this state the 
kinetic power of radio emitting outflows can significantly ex- 
ceed their radiative power both in AGNs and binaries (e.g. 
Owen et al. 2000; Gallo et al. 2003). The nucleus of M87 
can be regarded as the prototypical example of an AGN in 
the low accretion rate mode. In M87, the energetics derived 
from the observed jets and radio lobes suggests an outflow 
power on the order of 10 44 erg s _1 , requiring an accretion 
rate of order 3 x 10~ 4 M E dd onto a 3 x 1O 9 M black hole 
for an outflow efficiency (power/Mc 2 ) of 10%. On the other 
hand, the current X-ray luminosity of the nucleus is less than 
10 41 erg s _1 and even after bolometric correction it remains 
10-100 times lower than the outflow power. 

For binaries the F r /Fx correlation breaks down when 
higher X-ray fluxes are considered (Fender et al. 1999; Gallo 
et al. 2003). This break is associated with the increased mass 
accretion rate and a transition of the source to the so-called 
"high/soft" state. In this state the bulk of the observed X- 
ray emission comes from a standard, optically thick and ge- 
ometrically thin accretion disk (Shakura & Sunyaev 1973), 
and the spectra are reasonably well described by the mul- 
ticolor black body model. In this state the radio emission 
drops by a large factor (or completely disappears) suggest- 
ing quenching of the outflow. Once again one can identify 
a similar behavior for AGNs (Maccarone, Gallo & Fender 
2003), although the trend is less clear than for binaries. 

There is, therefore, good observational evidence that the 
energy release via accretion onto black holes in X-ray bina- 
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Figure 1. Sketch of black hole energy release as a function of 
mass accretion rate. At low accretion rates a significant fraction 
of accretion power goes into an outflow (thick solid line). Above 
a certain accretion rate, corresponding to the transition to the 
standard accretion disk mode, the outflow power decreases. The 
radiative power (thin solid line) is instead very low at low accre- 
tion rates and reaches a fixed fraction (~10%) of Mc 2 at accretion 
rates above 0.01-0.1 of Eddington. 



ries and AGNs depends on the accretion rate qualitatively 
as shown in Fig.l. At low rates there is a strong outflow 
whose power correlates with the accretion rate. The radia- 
tive efficiency at these rates is low, as e.g. in advection dom- 
inated flows (Esin, McClintock & Narayan 1997), so that 
the outflow dominates the output energetics. Numerous ob- 
servations of Galactic X-ray binaries have shown that at an 
accretion rate of order of ~ 10~ 2 — 10 _1 M E dd the sources 
make a spectral transition from the low/hard to high/soft 
state. We use this value as the characteristic accretion rate 
above which the outflow power decreases. In the particular 
case of GX339-4 the radio flux has dropped by a factor of 
more than 25 during the transition to the soft state (Fender 
et al. 1999), which can be viewed as evidence of a similarly 
strong decrease of the outflow power. In the absence of a jus- 
tifiable quantitative model of the outflow we adopt a simple 
prescription (shown in Fig.l) according to which the outflow 
power drops by large a factor when the accretion rate rises 
from ~ 10~ 2 M Edd to ~ M Edd . 

Despite the above justification of the picture presented 
in Fig.l, it may still be an oversimplification of the real 
situation. First, as is well known, ~ 10% of QSOs exhibit 
strong radio emission (e.g. Ivezic et al. 2002), indicating 
that powerful outflows may sometimes be associated with 
rapidly accreting SMBHs. What determines the radio loud- 
ness/quietness of QSOs remains an open question. It is pos- 
sible that apart from the instantaneous accretion rate, other 
properties such as black hole spin or the pre-history of the 
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accretion rate may play crucial roles. An indication that the 
latter may indeed be important comes from observations in 
transient X-ray binaries of a short, intense radio outburst 
during the transition from low/hard to high/soft state - 
in the so-called very high state - when the accretion rate 
approaches the Eddington limit (Fender, Belloni & Gallo, 
2004); such radio flares, however, have never been observed 
during the reverse transition (from high/soft to low/hard 
state). Importantly, it is not yet clear whether the observed 
radio outbursts are associated with a significant additional 
energy output or they mostly reflect a change in the emission 
properties of the previously (in the low/hard state) ejected 
material. 



3 BLACK HOLE FEEDBACK IN ELLIPTICAL 
GALAXIES 

We now discuss the conversion of energy released by accre- 
tion onto SMBHs into gas thermal energy in ellipticals. 

At low accretion rates the most important form is me- 
chanical energy input from outflows generated by the black 
hole. Examples of outflow interactions with the ISM are 
found in many elliptical galaxies (e.g. Finoguenov & Jones 
2001). If the outflow is halted well within the region occu- 
pied by the cooling gas, then the efficiency of converting the 
outflow energy into the gas heating can be high, up to 100% 
(e.g. Churazov et al. 2002). For simplicity we assume that 
the gas heating rate due to SMBH feedback is equal to the 
outflow power. This leads to a dependence of gas heating 
rate on the accretion rate (Fig.2) that mirrors the outflow 
power curve (Fig.l). 

Apart from the mechanical feedback, radiative gas heat- 
ing, via photoionization and Compton scattering of radia- 
tion from the SMBH, should be at work (Ciotti & Ostriker 
2001; Sazonov, Ostriker & Sunyaev 2004). Because of low 
efficiency of converting luminosity into heat the radiative 
feedback can be safely neglected at low accretion rates, but 
it can however be important at near-Eddington luminosi- 
ties (see below). It may therefore contribute significantly (or 
dominate) the high-M part of the heating curve in Fig.2 far 
from its maximum. This insures that total feedback power 
(mechanical plus radiative) does not go zero near Eddington 
rate even if the outflow quenches completely. 

It is convenient to characterize the gas heating power 
F through the efficiency <5 of converting accreted mass into 
heat: F = Sc 2 M. For subsequent analysis we assume that at 
very low accretion rates (outflow regime), the efficiency is 
high, So ~ const ~ 0.1, while close to the Eddington limit, 
the efficiency is much lower, Se ^ So- 

The schematic heating curve plotted in Fig.2 predicts 
the maximum feedback power at the level of 0.02 x LEdd, i.e. 
~ 2 10 45 erg/s for a 10 9 Mq black hole. Comparison with the 
most powerful outbursts identified so far in MS0735. 6+741 
(McNamara et al., 2005) and Hercules A (Nulsen et al., 
2005), both having an estimated outburst power at the level 
of ~ 2 10 46 erg/s, suggests that the maximum jet power can 
be an order of magnitude higher than adopted in Fig.2. Un- 
less these are very exceptional systems or their power is sig- 
nificantly overestimated, the efficiency of energy extraction 
by jets might need to be upscaled significantly from ~ 10% 
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Figure 2. Illustration of gas heating and cooling in elliptical 
galaxies. The thick solid line shows, as a function of the SMBH 
accretion rate, the heating rate due to outflow, which is com- 
plemented/dominated by radiative heating near the Eddington 
limit. Horizontal dashed lines show the gas cooling rate. The up- 
per cooling line represents a young galaxy in which a large amount 
of gas is present and/or the black hole is small. Feedback from 
the black hole is not able to compensate for gas cooling losses 
and the black hole is in the QSO stage with a near-critical ac- 
cretion rate, high radiative efficiency and weak feedback. As the 
black hole grows it moves down in this plot. The black solid dot 
marks the termination of this stage, when the black hole is first 
able to offset gas cooling, despite the low gas heating efficiency. 
The lower cooling line illustrates present day ellipticals: a stable 
solution exists at low accretion rates when mechanical feedback 
from the black hole compensates gas cooling losses. The radiative 
efficiency of accretion is very low and black hole growth is very 
slow. 



adopted here and/or the heating curve peaks at a higher 
accretion rate. 

Consider now a gas halo feeding a black hole at its cen- 
ter. If the black hole has a mass such that its maximum 
possible heating power exceeds the gas losses (so that the 
cooling line intersects the heating curve in Fig.2, e.g middle 
horizontal dashed line), the system can find itself in one of 
two equilibrium states (labeled A and B), in which cooling 
is balanced by heating. We argue below that only one of 
these states is stable. 

Consider for example that the mass accretion rate onto 
the black hole is described by the Bondi formula: 

-3/2 



M B ondi = 47tA(GM B h) c s p<x S 



(1) 



where A is a numerical coefficient that depends on the gas 
adiabatic index 7 (A = 0.25 for 7 = 5/3), c s = ^Jlj^- is the 

gas sound speed, p is the gas mass density, and s oc T e /p 7_1 
is the gas entropy index. Then it is the gas entropy near the 
black hole which regulates the amount of heat supplied by 



© 2001 RAS, MNRAS 000, 1-6 



4 Churazov, Sazonov, Sunyaev, Forman, Jones and Bohringer 



the SMBH into the gas. Consider the system near state B 
(Fig. 2), with a heating rate slightly exceeding cooling losses. 
The gas then heats up, its entropy increases, and it expands. 
The gas expansion will decrease the gas cooling losses. The 
Bondi accretion rate will also decrease. But according to 
our feedback prescription in the vicinity of B (see Fig. 2) the 
feedback power will instead increase. Thus the imbalance 
between heating and cooling will increase and the system 
will evolve away from state B. 

For the system in state A (Fig.2), the stability of the 
equilibrium depends on the specifics of energy deposition 
into the gas and on the gas distribution. We can verify the 
stability for a simple one-zone gas model. Consider a lump of 
gas with a fixed total mass M gas , kept at fixed pressure. The 
total cooling losses of the lump C oc n R 3 A(T) oc nA(T), 
where n is the gas density, R the size of the lump, and 
A(T) is the gas cooling function. If we neglect changes of the 
cooling function with temperature then C oc s~ 2 / 5 . This is 
a less steep function of entropy than the feedback power oc 
MBondi oc s~ 3//2 near state A and therefore the equilibrium 
will be stable. 



3.1 A Toy Model 

We now consider the implications of our model for the evo- 
lution of elliptical galaxies. We completely neglect stellar 
feedback, which can be especially important for low mass 
systems. Thus the discussion below mainly refers to massive 
systems with velocity dispersions larger than ~200 km/s. 

We follow the picture of a quasi-steady cooling flow in 
an isothermal halo (e.g. White & Frenk 1991). We assume 
that the potential of the galaxy is described by an isothermal 
sphere 4> = 2cr 2 lnr and the initial gas distribution follows the 
dark matter: 



Pgas 



fgas 



2tyGt 



(2) 



where a is the velocity dispersion, r is the distance from 
the center and f ga s is the gas mass fraction. The angular 
momentum of the gas is neglected. The initial mass of the 
black hole is assumed to be small enough that its feedback 
can be completely neglected (situation corresponding to the 
upper cooling line in Fig.2). The cooling radius r CO oi is set by 
equality of the cooling time t coo i = 3/2 l ( " n e .^ ) ) fcT « „^cr) 
to the age of the system t, where n e and Hi are the electron 
and ion densities: 



T"cool — 



fgasA(T)t 



1/2 



57 4 /2 A^ 2 kpc, 



(3) 



~ 23 erg s 1 cm 3 , 



where hereafter A(T) = A23 x 10" 

1T200 x 200 km/s, t = tg x 10 9 yr and the gas fraction is set to 
fgas = 0.17 (Cosmological baryon fraction). Inside the cool- 
ing radius the gas flows inward due to cooling losses, with 
the density of the gas increasing approximately as r~ 3 / 2 and 
the temperature being approximately equal to the virial one 
(T « 3 x 10 6 o- 2 00 K). At a radius r s <JC r coo i (at least for 
massive systems) the cooling time becomes comparable to 
the sound crossing time of the region, i coo i w r s /c s : 



fgasKT? 



9y3W 2 -ym^^aH 1 / 2 



1.7 a a 2 t g 1/2 A 3 fkpc. (4) 



In the vicinity of r a the gas cools rapidly and the tem- 
perature drops. We assume that if a simple cooling flow so- 
lution extends down to r s then cold gas will be present in 
large amounts in the system, stars will be actively forming, 
and the black hole will accrete at the Eddington rate, in- 
dependently of the details. The black hole mass then grows 
exponentially with the Salpeter e-folding time f s ?s4x 10 7 
yr: Mbh oc e'/* s . The mass cooling rate of the gas is instead 
a decreasing function of time: 



Mcool = 



2 /^VA(r) 1 / 2 

3? m p G :i / 2 t^ 2 



180 ^ootg 1/2 A 23 2 Ms/yr- (5) 



White & Frenk (1991) noted that the gas infall rate, calcu- 
lated as a time derivative of the gas mass within the virial 
radius (r vir = 2o km / s (l + z)~ 3 / 2 kpc), rather than the cool- 
ing rate may be a limiting factor. However for a redshift z of 
2-3 the infall rate is not much smaller than the cooling rate 
(see Fig.2 of White & Frenk 1991). Therefore one should 
not expect dramatic changes if the cooling rate is replaced 
by the minimum of the infall and cooling rates, especially if 
massive systems are considered. We also note that at 2 ~ 2- 
3 the cooling rate (Eq.5) is much larger than the Eddington 
accretion rate (for local Mbh I a relation). 

The gas cooling losses near r s are accordingly C w 



3/2fcT oc t 1 / 2 . Additional losses between r s and r coo i 
are of order of Cln(r coo i/r s ). With time cooling losses will 
decrease relative to the (growing) heating power of the black 
hole accreting at the Eddington rate. With our definition of 
the feedback amplitude, the black hole accreting at the Ed- 
dington rate will be able to offset gas cooling provided that: 



S E c 2 M Edd ^C^3/2kT 



Me, 



(imp 



(6) 



Substituting the expression for the Eddington accretion rate 
(oc Mbh) and solving for Mbh we find: 



M B 



4 x 10" 



4 ri/2 » m l 2 
V <y_ JgasHT) ' OT 

5e 2\/37r 3 / 2 cG 5 / 2 rn 2 i 1 /2 



(7) 
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for an accretion radiative efficiency r\ = 0.1. Once the 
black hole has grown to this size it will be able to pre- 
vent gas cooling, cold gas will disappear and star forma- 
tion will cease. The entropy of the gas increases dramat- 
ically and the system evolves toward a low accretion rate 
solution. In this regime the accretion rate required to pro- 
duce the same heating power as before is substantially sub- 
Eddington: M ~ |^M Ed d. 

We can infer from equation (7) that for a black hole 
mass satisfying the observed Mbh-o" relation (Tremaine et 
al. 2002), a balance between heating and cooling can be 
achieved for a quite low feedback efficiency (power/Mc 2 ) 
~10~°. The minimum possible feedback is provided by ra- 
diative heating, whose efficiency is (10~ -10~ )r (Sazonov 
et al. 2005), where r is the optical depth of the heated gas. 
Using equations (2) and (4), we can estimate the optical 
depth of the cooling flow external to r s as r w 0.004O-I00A23 1 . 
Furthemore, for relatively large galaxies (a > 200 km s _1 ), 
the cooling flow extends to a few kpc from the SMBH, 
where gas can be heated on the time scale of SMBH growth 
(Sazonov et al. 2005). Therefore, radiative feedback alone 
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appears sufficient to suppress the cooling flow once the black 
hole accreting at a near-critical rate reaches the observed 
Mbh-<t relation. Any additional channel of energy input into 
the cooling gas, e.g. via relativistic jets, will only faciliate the 
establishment of thermal equilibrium. 

The resulting system evolution is shown in Fig. 2 as the 
gray track. The black hole first accretes at the Eddington 
level (vertical track) and then switches (horizontal track) to 
the low accretion rate mode. The vertical track corresponds 
to a very high radiation efficiency and active star formation 
(cold gas is continuously deposited) , whereas in the low ac- 
cretion regime the radiative efficiency is low and little star 
formation takes place in the galaxy. 

Importantly, for the black hole mass for which the ther- 
mal equilibrium is reached, the maximum possible feedback 
power (maximum of the function F, see Fig. 2) significantly 
exceeds the cooling losses of the gas. This has two possible 
implications. First, the transition to the low accretion rate 
state might be followed by a very large energy input from 
the SMBH and ejection of gas from the system. Present day 
ellipticals are rather gas poor and thus may have suffered 
significant gas loss during the transition phase. On the other 
hand, the large maximum feedback power ensures that the 
system will be able to find a stable equilibrium at a low ac- 
cretion rate to support a large atmosphere of hot gas even 
without gas ejection. The latter case is probably relevant 
to cD galaxies in rich clusters. Since the gas heating effi- 
ciency in the low accretion rate state is high, the total mass 
accumulated by the SMBH over ~ 10 Gyr of passive evolu- 
tion will not be significant: Mbh ~ 2 x 10 7 Z/43 (0.1/So) Mq, 
where L43 is the gas cooling losses in units of 10 43 erg s _1 . 
Only for rich clusters can the black hole significantly increase 
its mass in this mode (e.g. Fabian et al. 2002) . 

In rich clusters the cooling losses might exceed the max- 
imum of the feedback curve. The central source can then en- 
ter a QSO mode again. The new episode of rapid accretion 
and black hole growth will continue until the SMBH grows 
enough to switch back to the low accretion rate mode. Such 
bursts of rapid accretion would be reminiscent of the models 
considered by Ciotti & Ostriker (2001). Similar episodes of 
QSO-type activity could be initiated by a merger with gas 
rich galaxies (as in simulations of Springel, Di Matteo, & 
Hernquist 2005) when large quantaties of cold/cool gas are 
suddenly dumped onto the existing black hole. 

The above model assumes that when the gas is able 
to cool the SMBH is accreting at the Eddington rate. 
The results would be different if the accretion rate is set 
by the Bondi rate from the hot phase of a multiphase 
medium (along the line of reasoning suggested by Nulsen 
and Fabian 2000). The thermal instability within r s , i.e. 
where the cooling time becomes shorter than the sound 
crossing time (e.g. Fall & Rees 1985) limits the density of 
the hot phase to maintain approximate equality of these 
time scales. As a result the Bondi accretion rate is limited 
to M « (t>j^rjcaTM Edd w <f> 0.3 (T/10 7 ) A^M E dd, where 
the dimensionless coefficient 4> of order unity depends on 
the detailed assumptions (see Nulsen & Fabian 2000). For 
the most massive systems, this maximum accretion rate is 
not far from the Eddington rate, while for low masses, the 
maximum rate may be several orders of magnitude below 
Eddington. 



4 SUMMARY 

Guided by the analogy of SMBHs and the stellar mass black 
holes in X-ray binaries, we postulated a non-monotonic de- 
pendence of the SMBH feedback power on the mass accre- 
tion rate associated with a transition from geometrically 
thick accretion flow to a geometrically thin disk. The gas 
heating rate is maximal at substantially sub-Eddington ac- 
cretion rates and drops strongly for near-Eddington rates. 
Such a scenario naturally leads to two stages in the evolu- 
tion of a massive elliptical galaxy: i) QSO-like nucleus and 
active star formation at early epochs and ii) weak nucleus 
and passive evolution at the present time. 

The main difference of the picture presented here with 
respect to previous work is that the phase of active black 
hole growth and star formation terminates when the SMBH 
produces enough power to substantially alter the thermal 
state of the infalling cooling gas, rather than to immedi- 
ately expell the bulk of the gas from the galaxy, which re- 
quires much more energy (see e.g Silk & Rees 1998; Wyithe 
& Loeb 2003). In our picture, heating the gas increases its 
entropy and also enhances the SMBH feedback power. The 
gas can (and probably will) be expelled but this may hap- 
pen over a long period of time after the end of the QSO 
stage. We showed that a Mbh oc ct 4 dependence with the 
observed normalization naturally arises for the simplest sce- 
nario of an isothermal gas halo if the SMBH feedback effi- 
ciency in the QSO phase is low (~ 10 -5 ) and not varying 
significantly with black hole mass. Furthermore, the postu- 
lated non-monotonic feedback curve results in "over-grown" 
black holes, i.e. SMBHs capable of supporting much more 
massive atmospheres of hot gas than those present at the 
end of the QSO stage and afterwards. 
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